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ABSTRACT 



In this paper we present the results of two detailed N-body simulations of the 
interaction of a sample of four massive globular clusters in the inner region of 
a triaxial galaxy. A full merging of the clusters takes place, leading to a slowly 
evolving cluster which is quite similar to observed Nuclear Clusters. Actually, 
both the density and the velocity dispersion profiles match qualitatively, and 
quantitatively after scaling, with observed features of many nucleated galaxies. 
In the case of dense initial clusters, the merger remnant shows a density profile 
more concentrated than that of the progenitors, with a central density higher 
than the sum of the central progenitors central densities. These findings support 
the idea that a massive Nuclear Cluster may have formed in early phases of the 
mother galaxy evolution and lead to to the formation of a nucleus, which, in 
many galaxies, has indeed a luminosity profile similar to that of an extended 
King model. A correlation with galactic nuclear activity is suggested. 

Subject headings: galaxies: nuclei, galaxies: star clusters, globular clusters: general, 
methods: n-body simulations 
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1. Introduction 



In the first paper of this series devoted to the study of the interaction of globular 
clusters (GCs) in triaxial galaxies, we analyzed the face-on collision between two GCs 
moving on quasi-radial orbits in the central galactic region, with the aim, also, to 
understand how effective is the tidal distortion. Our first finding was that the face-on 
collision is practically uneffective respect to the role of the external tidal field. Actually, the 
tidal erosion has been shown to be such to destroy loose GCs (initial King concentration 
parameter c < 1) in few passages through the galactic centre, while tight cluster (c > 1.6) 
keeps bound a substantial amount of their mass up to the complete orbital decaying. At 
this latter regard, another important result is that the orbital energy dissipation due to 
the tidal interaction is of the same order of that caused by dynamical friction. At the 
light of these results, and given that dynamica l friction was shown to be important in 



segregating massive GCs in triaxial potentials (IPesce. Capuzzo-Dolcetta &: Vietril Il992l 
Capuzzo-Dolcetta 1993 ; Capuzzo-Dolcetta fc Vicari 20051 ). a natural further step in the 



present program of investigation is to study the possible merging of a set of GCs decayed in 
the central galactic region to see whether a sort of 'super star cluster '(SSC) results from 
the merging, and what its morphological and dynamical characteristics could be. 

It is quite ascertained the existence of very bright (10 7 -r- 1O 8 L ) clusters of stars in the 
central region of galaxies across the Hubble sequence thanks to VLT an d HST observations 
(IWalcher et al J 120051 : Ivan der Marel et al.l 120071 : IWehner fc Harris! bood ) . On the basis of 
integrated colours and of the estimated M/L ratios, these clusters are usually thought to 
be young or, at least, to contain a significant population of young stars. In any case, that 
of the age is a controversal point, because, for instance, the nuclear star clusters (NCs) 
in M82 show evidence for mass segregation despite their spectra ar e well fitted by stella r 
population synthesis models with ages 10-50 Myr ( McCrady 2004 ). Boker et al.l ( 2004 ) 
fitted analytical models to Hubble Space Telescope images of 39 NCs in order to determine 
their effective radii after correction for the instrumental point-spread function. They 
compared the luminosities and sizes of NCs to those of other ellipsoidal stellar systems, 
in particular the Milky Way globular clusters, finding for NCs a narrow size distribution 
statistically indistinguishable from that of Galactic GCs, even though the NCs are, on 
average, 4 mag brighter than the old GCs. They discuss some possible interpretations of 
the similarity among NCs and Galactic GCs and, from a comparison of NCs luminosities 
with various properties of their host galaxies, find that more luminous galaxies harbor more 
luminous NCs. It remains unclear whether this correlation reflects the influence of galaxy 
size, mass, and/or star formation rate. 



However, iRossa et al.l (120061 ). by means of spectroscopic HST/STIS data, derived 
NCs ages for 19 galaxies and found that they range from 10 Myr to 10 Gyr, with a non 
negligible presence of old clusters. A c omparable result has been deduced also for 9 NCs 
of very late-type, bulge-less spirals, by IWalcher et al.l (120061 ) by means of a high-resolution 
spectroscopic survey with VLT/UVES. 
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Genzel et al. fl2003h . in their study of the stellar cuspy distribution around the 



Supermassive Black Hole (SBH) in the Galactic center, found that the K-band luminosity 
function of the local NC (within 9" of Sgr A*) resembles that of the large-scale Galactic 
bulge except for showing an excess of stars at K s < 14. It fits with population synthesis 
models of an old, metal-rich stellar population with a contribution from young, early, and 
late-type stars at the bright end. In the central arcsecond, they argued that a stellar merger 
model is the most appealing explanation. These stars may thus be "super blue-stragglers" , 
formed and "rejuvenated" through mergers of lower mass stars in the very dense (> 10 s M 
pc -3 ) environment of the cusp. 

Another intr i guing piece of the puzzle of the structure of SSCs is given by 



Baumgardt et al.l (120031 ). Through a comparison between the observational data on the 
kinematical structure of the very bright cluster Gl in M31 — obtained with the HST WFP 
Camera 2 and Space Telescope Imaging Spectrograph instruments — and their results of 
dynamical simulations carried out using the special purpose computer GRAPE-6, they 
obtained very good fits when starting simulations with initial conditions extracted from the 
end product of a previous simulation concerning a merging between two pre-existent star 
clusters. The merging explains observed features without the need to invoke the presence 
of an intermediate- mass black hole in the center of Gl. 

Not many simulations have been presented in the literatur e that study the possible 



form ation scenarios for SSCs. Among these, we remind those by lFellhauer &: Kroupal (12002 



20051 ) finding that star clusters aggregates, like the ones found in the Antennae or Stephan 
Quintet, are very likely the merger progenitors of SSCs. Interestingly enough, these authors 
also find that the resulting SSC is a stable and bound object, whose density profile is well 
fitted by a King profile, even if the mass loss of the merger product occurs through every 
perigalacticon passage. 



2. Globular cluster merging and galactic nuclei 



The increasing amount of data about massive clusters in galaxies, together with that 



of GC systems in galaxies, especial ly of early type but also for spirals (jHarria Il986l . 11991 



Ashman fc ZepJ Il998l ; lHarris 



2001, etc.) makes interestin g the inv estigation, started in 



the first paper of this series (jMiocchi et al.ll2006l hereafter iPaper J ), of the fate of massive 
GCs moving i n the par ent galaxy field, subjected to dynamical friction braking and tidal 
interaction. In IPaper II we studied the dynamical evolution of two GCs with mass ~ 10 7 M 
moving on quasi radial orbits in a triaxial galactic potential, following eight passages across 
the galactic center. The scope of that paper was to investigate the chances of survival 
of GCs to possible tidal disruption induced both by the external field and by mutual 
cluster-cluster interaction; with this aim, we maximized these effects by giving the clusters 
initial conditions corresponding to quasi-radial orbits. 
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In this Paper II, we have another aim, that is to study whether and how the merging of 
various massive GCs decayed by dynamical friction in the inner galactic region may occur. 
The main questions to answer are: (i) given some (realistic) initial conditions for a set of 
GCs which experienced a significant orbital decay, are they undergoing a full merge? (ii) if 
so, what is the time needed? (iii) what is the final structure of the merged NC? (iv) does it 
attain a quasi-steady state? 

The answers to these questions are of overwhelming importance to give substance 
to the interpretation of the formation of early type galaxy nucl ei via merging of 



extensively studied by 


CaDuzzo-Dolcetta 


Capuzzo-Dolcetta & Vicari ( 


2005 


)■ 



( 1975 1) and subsequentl y 



This scenario of galactic nuclei formation has been raised again recently, due to 
the independent and almost contemporaneous findings that NC masses obey to sim ilar 



scaling relationships with host galaxy properties as SBHs do (see iRossa et al. 



2006 for 



2006 for 



spirals; IWehner fc Harris! 120061 fo r dwarf elliptical [ dE] galaxies and ICote et al. 
elliptical galaxies). In particular, ICote et al.l (120061 ) give evidence, in early-type galaxies 



in the Virgo cluster, of NC luminosity distributions that a re much better fit ted by an 
extended profile (King's) rather than by a point source (s ee lBoker et al.ll2002l for a similar 
finding in late-type spirals and iGraham &: Guzman! 120031 for dE galaxies in the Coma 
Cluster). The half mass radii of nuclei (r/J are in the range 2 < r^(pc) < 62, with 
(r/j) = 4.2 pc and correlate with the nucleus luminosity: oc L° n 5±om . The mean of 
the frequency function for the nucleus-to-galaxy luminosity ratio in nucleated galaxies, 
log 7/ = —2.49 ± 0.09 is indistinguishable from that of the Super Massive BH-to-bulge mass 
ratio, \o£{Mbh I M nai) = —2.61 ± 0.07, calculated in 23 early-type galaxies. On these bases, 
Cote et al.l (120061 ) argue that resolved stellar nuclei are the low-mass counterparts of nuclei 
hosting SBHs detected in the bright galaxies. If this view is correct, then one should think 
in terms of central massive objects, either SBHs or SSCs, that accompany the formation 
and/or early evolution of almost all early-type galaxies. 

It is clear that these characteristics of galactic nuclei well fit into a scenario of multiple 
GC merging in the inner galactic regions (that we call "dissipationless" scenario) as 
alternative to the "di ssipational" scenario , this latter being based mainly on s peculative 
hypotheses (see, e.g.. Ivan den Berghlll986l: ISilk et al.l Il987l : Babul fc Reeslll992h supported 
by some quantitative results (IMihos fc Hernquistl Il994j ). The proof of the validity of 
the "dissipationless" hypothesis requires, at first, a detailed iV-body simulation of the 
interaction of stellar clusters in the inner region of their parent galaxy, taking into account 
both the cluster-cluster and the cluster-galaxy dynamical interaction. This latter includes 
tidal distortion, acting on the cluster internal motion, and dynamical friction, acting on the 
cluster orbital motion. 



In this context, we cite the encouraging results of the numerical simulations made 
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by lOh fc Linl (120001 ) who revisited the hypothesis of dE nuclei formation through the 
orbital decay of GCs and suggested that this occurs mainly in galaxies with a relatively 
weak extragalactic tidal perturbation, leading to the formation of compact nuclei within 
a Hubble timescale. They also show that the central galactic field does not destroy the 
integrity of the clusters and facilitates the merging to occur. Moreover, they find that 
the observed central structures of some nucleated Virgo Cluster dEs are well reproduced 
by superimpo sing a sma ll number of globular clusters to the galactic stellar distribution. 



Nevertheless, lOh fc Lin! used a very small number of particles (= 500 in each cluster) to 



simulate the final merger stage, during which they neglected the role of dinamical friction, 
as well. Thus, their results could be affected by spurious collisional effects, even though a 
rather large smoothing radius was adopted in the inter-particle interactions. 



More recently, iBekki et al.l (120041 ) show how the formation of NCs via multiple merging 
of GCs, leads to systems with global scaling relations which are compatible with those 
observed for galactic nuclei. Unfortunately, as mentioned before, these results cannot be 
considered conclusive because in these authors' simulations the, important, role of the 
external field was neglected. 



3. The simulations 

We consider GCs as iV-body systems moving within a triaxial galaxy represented by an 
analytical potential, subjected also to the deceleration due to dynamical friction (hereafter 
df). To keep a high level of resolution in the simulations (i.e., a large number of particles 
to represent each GC without exceeding computational capabilities, we decided to limit our 
study to the interaction of a limited number, four, of GCs. 

We studied the merging process of these GCs, as if they were already decayed to the 
inner region of the galaxy, in two separate simulations (cases 1 and 2). A quick orbital 
decay induced by df is due to initial large values for the total mass of the clusters. These 
large masses are compatible with those of many young GCs actuall y obsery ed in various 



galaxies (e.g.. iKissler-Patig et al.ll2006l . and references therein), see iPaper II . for a deeper 
discussion. The 4 clusters are labeled (al), (bl), (cl), (dl) in the case 1 and (a2), (b2), 
(c2), (d2) in the case 2. In the simulation of case 2, the 4 GCs have the same initial orbital 
conditions but their scaling parameters are such to give them less dense and more extended 
initial configurations (see also Tableland Sect J3.2I) . 



Unless otherwise specified, the same units of measure used in IPaper t are adopted in 
this paper: lengths, masses and time are measured, respectively, in unit of the galactic 
core radius rb, of the galactic core mass Mb and of the galactic core-crossing time 
tb = [rl/ '(GMb)] 1 / 2 . Note that the actual physical values for these quantities are unrelevant 
because the results of simulations can be scaled as long as the 2-body stellar collisions 
are negligible; otherwise, the dynamics would obviously depend also on the mass of the 
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individual star in the clusters. 



As regards the computational techniques, we adopted an iV-body representation for 
the stars in the clusters, simulating their dynamics by means of the parallel 'tree ATP' co de 
(jMiocchi fc Capuzzo-Dolcettall2002l ). whose main features were also described in lPaper I . 



3.1. The galactic model 



We adopt the same galactic model u sed in the simulations of iPaper 1 i.e. the 



self-consistent triaxial model described in Ide Zeeuw &: Merrittl (119831 ) . It corresponds 



to a non-rotating ellipsoidal and triaxial distribution of matter oriented in such a way 
to have the longest and the shortest principal axis aligned with the coordinate axes x 
and z, respectively. The axial ratios are 2 : 1.25 : 1, leading to a projected profile in 



agreement with that observ ed in the spheroids of some spirals ( see, e.g., iBertola et al.lll 991 



2001 



Matthews fc de Griia I2004T) and in elliptical galaxies (see, e.g., IWagnerlll988l ; iDavies et al 



Bertola et al. 


1991 


1988 




Davies et al. 



Statler et al.l I2004J ) . The resulting potential can be expressed as the sum of a 



spherically simmetric density following the modified Hubble's law 

-3/2 

f r \ ~ 

Pb(r) = pbo 

■ r b 



with scale parameters (core r adius) and Pbo = My,l r\, plus other two non-spherical terms 
that give the triaxial behaviour (jSchwarzschildl Il979l ) . The mass parameter Mb is the mass 
enclosed in a radius slightly smaller than r&, in fact Mb ~ 0.45M(rfe). The mass in the 
generic sphere of radius r is contributed only by the spherically simmetric part, Pb{r), of 
total density, giving: 

-2 \ 1/21 



M(r) = 47rM b < In 



r 

— + 



r 

^7 + 1 



+ 1 



-1/2 - 



(2) 



that gives an infinite total mass. However, there is no need of a cutoff in the model, because 
only the gravitational potential (which is not divergent) is important for the purpose of our 
simulations. 



As in IPaper II . the reference frame has the origin at the galactic center and the x and z 
axes are, respectively, along the maximum and minimum axis of the triaxial ellipsoid. The 
centers-of-mass (CM) of the clusters were initially located well within the galactic core (see 
Table E]). 



With regard to the df, we used the generalization to the triaxial case (see iPesce et al. 



19921 ). of the Chandrasekhar formula (jChandrasekharl Il943l ). with a self-consistent 
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evaluation of the velo city disp ersion tensor, taking also into account that the GCs are 
extended objects. See IPaper I for more details at this regard. 



3.2. The cluster model 



Our numerical study involved two sets of four clusters evolved in two separa te 



simulations. The GCs initial internal distribution was sampled from a stellar iKingi (119661 ) 
isotropic model with total mass M, velocity parametei0 <tk, limiting and "King" radius rt 
and r c , respectively; c = \og(r t /r c ) is the concentration parameter, t ch = [r^/ (GM)] 1 / 2 the 
half-mass crossing time and po is the central density. The 'limiting radius' is the radius 
at which the King distrib ution function drops to zero, to reproduce the presence of the 
external field (lKingill966l ). The initial values of these parameters are listed in Table HJ 



Each cluster was represented with N = 2.5 x 10 5 'particles', whose individual masses 
were assigned according to a Salpeter's mass distribution (dN oc m~ 2 - 35 dm) cutted at m min 
and m max = 100m min , with average (m) = 3.1m min , where m min is adjusted to give the 
desired cluster total mass (see Tabled]). When simulating massive GCs, the assumption 
N = 2.5 x 10 5 gives a high value of m m i n that shifts the chosen mass distribution in form 
of Salpeter law towards large masses so t o lose it s representativity for a real GC stellar 
mass spectrum. However, as discussed in IPaper J . we use unequal mass particles to verify 
that no collisional relaxation takes place during the simulation, as demanded by the 2-body 
relaxation time of the clusters which is in any case longer than the time to merge. 

Actually, even for the most compact case considered here (model (al) see Table [Q, the 
half-mass relaxation time of the simulated cluster is larger than the merging time t m ~ 18 
(see Sect. H|) after which the clusters lose their individuality and evolve as a larger N system 
(see discussion in SectJSJ) The real cluster (with a larger number of stars) relaxation time is 
obviously longer than that of the simulated (sampled) cluster. 

Were the GC s massiv e enough to decay rapidly into the galactic core, it can be safely 
assumed (see also lPaper J . sect. 2) that their age at the beginning of the simulation is less 
than their internal 2-body relaxation time, so that mass segregation has not occurred and, 
thus, it was not included in the initial stellar distribution. 

In the following, we often refer to the system center- of- density (CD) - i.e. the average 
particle position weight ed with the loca l density instead of the mass - as defined in 
Casertano fc Hutl (119851 ) . As in IPaper t . in most cases we took the CD as the origin of the 
best suitable reference frame for the study of the internal system properties. 



1 Which is such that, at the center , a 2 = (v 2 ) ~ 3cry, the ap proximation being as better 



as the model concentration is higher (IBinney fc Tremaindll987l sect. 4.4c). 
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4. Results 

In spite that the true conserved quantity, in absence of any external dissipation, is 
the center-of-mass orbital energy, it is preferable to refer to the CD orbital energy instead, 
because the CD identifies better the actual GC location when the outer part of the cluster 
is highly distorted and dispersed. Defining E orb as the orbital energy (per unit mass) of 
the cluster CD, its dissipation due to both the df braking effect and to the tid al interaction 
with the environment, can be quantified by the following parameter ( Paper j ): 



_ E orb (t) - fr 

Mt) ~ £W0)-V (3) 

where \l/o = —4:iiGMb/rb is the central galactic potential well. In Fig. [3j £ rb(t) is plotted 
for all the clusters in the two cases studied. 

The merging occurs rather quickly, as it is indicated by the snapshots of the system 
as projected onto one of the coordinates plane (see Fig.s [1] and [2] for the case 1) and 
clearly seen from the time behaviour of the Lagrangian radii of the whole system (Fig. HJ). 
Lagrangian radii attain steady values at a "merging" time of about t m ~ 18, when also the 
innermost radius shown stabilizes; note, that the inner is the Lagrangian radius, the later is 
the time at which the stable state is achieved. A similar merging time can be deduced also 
by the time behaviour of the distances of the clusters CD to the galactic center (Fig. |3J). 

Of course, the energy dissipation rate is closely related to the prescrip tion adopted for 
the evaluation and treatment of df. For instance, our evaluation ( Paper 1 Appendix A) 



could lead, in principle, to an overestimate of the effect, because we consider the cluster as 
a constant mass pointH Obviously, df would be accounted for accurately when adopting a 
full iV-body representation of the galactic region in which GCs move; for this reason we are 
running a self-consistent simulation that we will present in a forthcoming paper. However, 
preliminary results show a clearly shorter orbital decaying time for the same clusters of 
case 1, but moving in a self-consistent representation of the same galaxy sampled by 5 x 10 5 
"particles" . This is shown by the time behaviour of the distance of the clusters CD from 
the galactic center depicted in Fig® from which a merging time t m ~ 13 is deduced, i.e. a 
factor 1.4 smaller than in the simulations with the analytic df form. This finding agrees 
with some recent results that show how stars stripped from the cluster by the field, but 
still close enough to th e system, continue to contribute to mass of the decelerating system 



( iFellhauer fc Lin! 120071 ). Moreover, recent fully self-consistent iV-body simulations confirm 



that the real df effect can be even stronger that that estimated b y the usual Chan drasekhar 



formula (probably for the further friction due to tidal effects, see lFujii et al.ll2007l ). 



2 To partially overcome this problem we conservatively used a halved cluster mass as mass 
parameter in the df formula. 
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4.1. The NC morphology 

To study the NC morphology, we evaluated the eigenvalues {Ii,l2,h} and the 
eigenvectors {ei,e 2 ,e 3 } of the inertia tensor 

1 N 1 

^ = ~N '^2( X i ~ X i)( X j ~ X °j) ~ Tp X 
k=l 

N N 
k=l k=l 

with (x\, x%, £§) being the position vector of the fc-th particle, and (x^x^x®) the CD 
position. In Fig. [7J the square root of the eigenvalues are plotted as a function of time. 
They are evaluated considering, in Eq. (TJ]), only particles closer than r h to the system CD, 
where r h is the half-mass radius of the NC at the end of the simulations. Since \JTi is equal 
to the length of the i-th axes of the ellipsoid fitted to the matter distribution, then their 
ratios give a direct measure of the triaxiality of the system; e.g. if 1\ — l-i — h then the 
system is spherical. 

It can be seen that, after the initial merging phase, the NC as a whole attains a stable 
configuration. Only in the case 2 (Fig. [7J upper panel) the system shape is moderately 
axisymmetric (oblate) and, by examining the orientation of the eigenvectors, it was found to 
have principal axes nearly aligned to those of the galactic potential, with axial ratios \fl{ : 
\JT~2 '■ \Zh — 1.4 : 1.4 : 1, and with e3 nearly parallel to the z-axis. This alignment suggests 
that the final morphology is influenced, to some extent, by the galactic morphology. Thus, 
it is not surprising that in the NC 2 outskirts the stars tend, eventually, to move on orbits 
compatible with the underlying galactic triaxial configuration, as we will see from the 
velocity anisotropy (Sect. 14.1.21) . Nevertheless, the self-gravity of the simulated system is 
important enough to make the shape of the NC 2 closer to a spherical configuration than the 
galaxy. Accordingly, the configuration of the NCi is almost spherical due to its stronger 
self-gravity (Fig. [7J lower panel). 



4-1.1. Density profile 

The radial density profile of the NCs, which we have seen to be quasi-spherical at the 
end of the merging, is plotted in Fig. M, where, for comparison, the profile corresponding 
to the spatial superposition of the 4 progenitor clusters, in their initial configurations, is 
reported, too. This gives an immediate 'visual' indication of the efficiency of the merging 
process; the merging remnant density profile is even more concentrated than that given by 
the mere sum of the initial profiles of the four progenitors. Actually, the merger profile 
has a smaller core (showing a central density a factor 1.25 higher) and a more extended 
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envelope. In fact, it is remarkably well fitted by a high-concentration (c = 2.2) King profile 
(see Fig. [8]). This result — w hich ag r ees wi th what found in the different context of merging 



galaxies since the pioneering IWhitd (119781 ) simulations — suggests that a sort of violent 
relaxation took place during the merging. This relaxation is likely due to the rapidly 
varying potential acting on each star of the clusters. 

In the initial cluster models a stellar mass spectrum was included without mass 
segregation, in order to investigate possible relaxation phenomena due to either spurious 
collisional effects or induced by the external tidal field. The analysis of the behaviour of 
the average mass as a function of the distance from the center in the final NC configuration 
indicates no evidence of mass segregation, thus enforcing the above-mentioned hypothesis 
about the violent relaxation as a cause of the increased concentration of the merging 
product. 



4-1-2. Velocity distribution 



To compare our simulation results with observational data giving the radial profile 
of the line-of-sight velocity dispersion, we must evaluate the "total" velocity dispersion, 
°~toti by properly accounting for the contribution of the underlying galaxy, this latter being 
given d irectly b y the self-consistent model (according to the fitting formulae [A3] and 
[A4] in iPaper J ). Fig M shows clearly a behaviour of a to t decreasing towards the galactic 
center, where the surface density of the NC dominates and thus o~ tot lowers to the velocity 
dispersion of the NC, which is 'colder' than the host galactic core. The behaviour shown in 
Figj9] may seem peculiar for a self-gravitating system, which normally keeps its equilibrium 
through a negative (or null, if isothermal) gradient of a 2 . This apparent peculiarity is due 
to that the central "observed" velocity dispersion is dominated by the NC which is not a 
self-gravi tating system, ind eed. This feature was actually found by a Keck II spectroscopic 
analysis (jGeha et al.ll2002l . see their Fig. 5) for most of the Virgo cluster nucleated dEs 
of the sample. The same atot radial behaviour is also show n by the solution of the Jeans 
equations for a sample of NCs observed in late- type spi rals (|Walcher et al.ll2005l ) and is also 
deduced on the basis of direct iV-body simulations by lOh Sz Lin! (120001 ) concerning GCs 
merging at the center of dEs. This latter simulation confirms the occurrence of a merging 
event within the central potential well of the host galaxy, even if care is due to the low 
number of particles that forces them to adopt a rather large gravitational smoothing radius, 
in order to suppress fictitious collisional effects. 

As regards the velocity anisotropy, we studied the anisotropy parameters (3^ = 
1 — a 2 e /a 2 of the stellar velocity distribution in the NCs, as a function of the distance from 
the system CD, where and 9 denote the azimuthal or the polar angle, respectively. These 
parameters are compared with those of the host galaxy. 



We have found that NCi has a nearly isotropic velocity distribution < 0.2), while 
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the influence of the external galactic potential induces a < (3q < (i.e. a^> Oq > of) 
in the outer region (r > 0.5r&) of the less compact NC2, which means that the velocity 
distribution in the outskirts is pre ferentially t a ngent ial and parallel to the xy plane, as 



in the galactic model. Noticeably, iGeha et al.l (120021 ) found that the a to t decrease in the 



center of their observed nucleated dE can be reproduced, in most cases, by a model having 
a nearly isotropic or moderately tangential velocity distribution. 

Another relevant parameter is the central value of p/cr 3 , which is proportional to the 
density in the phase-space. Its value (see Table W) in the merged systems is smaller than 
that of the four individual progenitor clusters; this corresponds to the, expected, rarefaction 
in the phase-space of the inner part of the system after the mutual interaction of the 
merging clusters and, mainly, that with the external field. This time-dependent interaction 
leads to a certain diffusion in phase space even if the four GCs constitute a collisionless 
systems. Note how the rise of the central space density respect to the linear summation of 
individual GC densities in the case of NC2 is strongly balanced by expansion in velocity 
space (a 3 of NC2 is ~ 5 times higher than the average velocity dispersion of the progenitors) 
so to decrease for a factor 0.4 the phase-space density. 

As regards the virial equilibrium of NCs, it is worth noting that the King model fitting 
of NCi (FigJHJ) gives a velocity parameter cr K = {A^GpQ) l ' 2 r c /?> — 0.36, significantly lower 
than the true central value of the line-of-sight velocity dispersion ~ 0.52. This discrepancy 
is due to the non negligible influence of the external galactic potential, which kinematically 
"heats up" the NC. This has an important implication on the reliability of the NC mass 
determination by means of the virial relation M = rjapTh/G, with a p the projected velocity 
dispersion. The true mass value for our NCs is obtained with a virial parameter rj ~ 2 that 
is about 5 times lower than the usually adopted values based on the unrealistic assumption 
of isolated King mo dels. This finding is in remarkable agreement with the results of 



Walcher et al.l (120051 ) on the dynamical mass estimate of NCs observed in bulge-less spiral 



galaxies. 



5. Are resolved galactic nuclei compatible with a merging origin? 



An important result of this work is that a quasi-stationary NC may form as merger 
product of orbitally decayed massive clusters. The projected density profile in the 2 cases 



Boker et al. 


( 


2002); 


Geha et al. 


( 


2002): 



Cote et al.l (120061 ) is clear; of course, the ratio between the central surface total density and 
that of the galaxy central plateau depends on the number of merged GCs. This ratio is 
gj 18 for the four merging GCs of case 1 and ~ 1.5 for case 2. An analysis of Fig. 4 in 



Cote et al.l (120061 ) indicates ratios of the central luminosity and the inward extrapolated 
galactic luminosity in the interval 1 < L /L ga i < 100. In the assumption that the density 



13 



contrast scales linearly with the number of GCs merged (N m ) the range of the observed 
ratio corresponds to the interval < N m < 22 for case 1 and < N m < 267 for the 
case 2. Of course, the linearity is not guaranteed; actually, the modes of the merging of 
a more abundant sample of GCs deserves a careful, forthcoming, study. Moreover, a fully 
self-consistent study of the feedback among the merging GCs and the surrounding stellar 
field is needed to confirm the modes and time of orbital decay and merging. This study 
will constitute a check of the reliab ility of the dynamical friction treatment done here in 
the frame of the iPesce et al.l (119921 ) scheme, which, if proved, will result on the possibility 
of saving a huge amount of CPU time in future simulations. 



6. The merger remnant stationariety 

In both the cases studied here the merging occurs at ~ 20 t b and the resulting NC 
keeps its characteristics almost unchanged over the remaining simulated time (other 25 
t b ). This does not imply, necessarily, the secular stationariety of the merged system, which 
should be investigated extending sufficiently in time the iV-body simulations. Anyway, 
simple considerations on the relaxation times convince ourselves that the lifetime of the 
me rged system should be quite long. Actually, the half mass relaxation time of the cluster 



is f 1Spitzerlll987f ) 



N ( M, 



t:h = linN ) {V b ) * < 5 > 

which gives t rh = 1.1 x 10%, for the values N = 10 6 , r h /r b = 0.12 and M/M b = 6.4 x 10~ 2 
corresponding to the four merged clusters of case 1, and t r h = 2.6 x 10 3 t b , for the four 
merged clusters of case 2, where the only different value is that of the ratio rh/r b which 
results 0.21 in this case. In physical units, this means t r h = 5.8 Gyr and t r h = 13.5 Gyr 
in the 2 cases, respectively, assuming r b = 200 pc and M b = 10 9 M Q . Consequently, the 
merged system should stay almost stationary for a time much longer than that of our 
present simulations and, more relevantly, of the order of the age of galaxies, especially when 
thinking that the 'true 'number of stars in the system is from 64 to 6.4 x 10 3 times the one 
used in the simulations for M b in the range 10 9 -10 n M©, so to scale all the relaxation times 
estimated above for a factor in the range 50-3, 900. If other clusters merge to the center, N 
increases, as likely does, increasing furtherly t r f t . In conclusion, a collisional evolution for 
the nuclear cluster seems unlikely, but larger scale instabilities cannot be ruled out. So, the 
dynamical evolution of the merger remnant deserves a more extended in time simulation. 

In the, likely, hipothesis that the NC keeps its stability during subsequent merging 
events, the central density po growth shown by FigJH] can be reasonably scaled with the 
number N m of merging clusters, giving p ~ 575N m (case 1) and p ~ 3.7N m (case 2) 
in units of the central galactic density, p b0 . With M b = 3 x 10 9 M Q and r b = 200 pc as 
typical values of normal elliptical galaxies, one has p bQ = 375 M pc -3 ; consequently, for 



1/2 / x 3/2 
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the NC to reach the high density of the order of 10 7 M pc~ 3 , which is typical of the 
environment of massive accreting objects, N m ~ 46 (case 1) and N m ~ 7150 (case 2) are 
needed. Even with all the necessary cautions suggested by the various assumptions done, 
these numbers (particularly for case 1) are indicating that a NC formed by merging of 
individual sub-structures may actually have reached central densities large enough and 
in a time short enough to provide mass for an efficient accret ion of a possible massive 



galactic object, as suggested in various previous studies (e.g.. ICapuzzo-Dolcettal Il993 



Capuzzo-Dolcetta fc Tesserilll999l ; ICapuzzo-Dolcetta fc Vicaril 120051 ) . In light of this, it is 



clear the importance of further, more detailed and complete, study of the actual modes of 
formation and evolution of a NC via the merging process (Capuzzo-Dolcetta Sz Miocchi, in 
preparation). 



7. Conclusions 



In this paper we studied the modes of interaction of a sample of few (4) globular clusters 
whose orbital motion is limited to the core of the galactic triaxial field, because they have 
experienced a significant and rapid dissipation of the orbital energy by dynamical friction 
acting on their large initial masses. We studied 2 sets of 4 GCs of different initial conditions, 
characterised by a higher (case 1) and lower (case 2) central density. Dynamical friction was 
properly taken into accoun t by the triaxial generalization of the classical Chandrasekhar 
formula (IPesce et al.lll992l ) and it is the main cause of orbital energy dissipation up to the 
beginning of the merging process, when the individual cluster sizes are comparable with 
the orbital size. Since this stage on, the main responsible of the residual orbital energy loss 
is the tidal torque. The merger is completed in ~ 18 galactic core crossing times, i.e. in 
a time much shorter than the Hubble time and, in any case, small compared to the total 
orbital decaying time. 

The merged system keeps some of the characteristics of the preexisting objects, 
attaining a relaxed structure which has, in the case 2, a mildly triaxial shape (inherited 
by the environment) and a halo which is diffusing in the external field assuming its phase 
space properties. The case 1 results into a merger configuration which conserves better the 
spherical simmetry of the 'progenitors' and whose radial density profile is more concentrated 
to the center than that expected on the basis of the mere space superposition of the 4 
progenitors. This is, likely, the consequence of some violent relaxation of the merged cluster 
to the state of a quasi-stationary Super Star Cluster, whose stellar density distribution 
maintains a King shape and is very similar to Nuclear Clusters observed at the center of 
many galaxies. Even in the limits of our simulations (both on the number of merging objects 
and on the time extension of the integration), we may infer from a comparison of our results 
with the characteristics (surfac e brightness profile , integrated light and velocity dispersio n) 



ff - 

of the nuclei of many galaxies flBoker et all 120021 ; iGeha et al.l 120021 ; iWalcher et all 120051 



pc 
2( 



Cote et alJl2006h that such nuclei may have actually formed by the merging of few tens 
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(when merging progenitors are quite compact) up to few hundreds (when progenitors are 
looser) massive GCs decayed orbitally into the inner galactic regions. 

As remarkable result, the radial profile of the line-of-sight velocity dispersion of the 
merger remnant shows the same minimum at the galac tic center as actu ally found by 
observa tions of a sample of Virgo cluster nucleated dEs (IGeha et al.ll2002l ) and of late- type 
spirals ( Walcher et al.ll2005l ). Finally, an important theoretical result is that high densities 
(> 10 7 M pc~ 3 ), critical to allow a significant accretion onto a massive object, may be 
reached by mean of a merging of less than 50 clusters of the type considered in case 1. This 
may have important implications on the role of NCs in supplying the innermost galactic 
activity. 
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Fig. 1. — Snapshots of the simulation in the case 1 (projection on the x-y plane) from t — 
to t = 20.3. 
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Fig. 2.— Continuation of Fig. □ from t = 21.6 to t = 42. 
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Fig. 3. — Time behaviour of the fraction of the orbital energy kept by the clusters CD. Solid 
line: clusters (al), (a2); dotted: (bl), (b2); short-dashed: (cl), (c2); long-dashed: (dl), (d2). 
Left panel: NCi progenitors; right panel: NC 2 progenitors. 
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Fig. 4. — Evolution of the Lagrangian radii evaluated with respect to the global CD. They 
refer to 10%, 30%, 50% and 90% of the total mass of the whole system. Left panel: case 1; 
right panel: case 2. 
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Fig. 5. — Time behaviour of the distance of the clusters CD to the galactic center. Left 
panel: NCi progenitors; right panel: NC2 progenitors. Symbols are as in Fig. [3] 




Fig. 6. — As in Fig. [5]for the clusters of case 1, but moving in a particle self-consistent model 
of the galaxy (from Capuzzo-Dolcetta & Miocchi, in preparation). 
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Fig. 7. — Square root of the eigenvalues of the inertia tensor (solid line: Ji, dotted: I2, 
dashed: J 3 ) evaluated on the part of the system lying in the sphere centered at the CD of 
the whole system and with radius equal to the final NC's half-mass radius. Lower panel: 
case 1; Upper panel: case 2. 
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Fig. 8. — Radial volume density profile of the NC at the end of the simulation (solid line), 
compared with that corresponding to the sum of the initial density distributions of the 4 
clusters (dotted line). Left panel: case 1, right panel: case 2. The dashed line represents the 
King best-fit to the NC profile. The structural parameter of the King profile are reported. 
The distance r is to the galactic center. The vertical dotted line marks the NC half-mass 
radius. 
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Fig. 9. — Total 3-d velocity dispersion profile. Left panel: case 1; right panel: case 2. 
Dashed curve is the velocity dispersion of the galactic model only. The total a is computed 
by summing the galactic model velocity dispersion plus that of the NC, weighted by the 
surface density. 
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Fig. 10. — Projected surface density profile for the last NC configuration, summed to the 
galaxy background profile (long dashed line) so to give the total surface density (solid line). 
Left panel: case 1, right panel: case 2. 
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Table 1. Clusters initial parameters. 



cluster model 


M 






c 






r h 




Po 


tch 




o- K 




(al) 


1.5 x 10" 


-2 


0.16 


1.2 


1.1 x 10- 


-2 


2.3 x 10- 


-2 


770 


2.6 x 10" 


-2 


0.36 


3.2 x 10 3 


(bl) 


1.7 x 10" 


-2 


0.16 


1.0 


1.5 x 10- 


-2 


2.4 x 10" 


-2 


410 


2.8 x 10" 


-2 


0.36 


1.7 x 10 3 


(cl) 


1.8 x 10~ 


-2 


0.14 


0.99 


1.4 x 10- 


-2 


2.2 x 10- 


-2 


610 


2.3 x 10" 


-2 


0.41 


1.7 x 10 3 


(dl) 


1.4 x 10- 


-2 


0.14 


0.89 


1.8 x 10- 


-2 


2.4 x 10- 


-2 


250 


3.3 x 10" 


-2 


0.34 


1.2 x 10 3 


(a2) 


1.5 x 10- 


-2 


0.77 


1.2 


5.3 x 10- 


-2 


0.11 




6.1 


0.30 




0.16 


290 


(b2) 


1.7 x 10- 


-2 


0.78 


1.0 


7.3 x 10- 


-2 


0.12 




3.3 


0.31 




0.16 


150 


(c2) 


1.8 x 10- 


-2 


0.68 


0.99 


6.9 x 10- 


-2 


0.11 




4.9 


0.26 




0.18 


160 


(d2) 


1.4 x 10- 


-2 


0.72 


0.89 


9.3 x 10- 


-2 


0.12 




2.0 


0.37 




0.15 


110 



Note. — Parameters list for the initial cluster models, expressed in galactic units. Reported are: the 
GC mass (M), the limiting radius (r t ), the King concentration coefficient (c), the King radius (r c ), the 
half-mass radius (r h ), the central density (p ), the half- mass crossing time [t ch = [r\/ (GM)] 1 / 2 ], the King 
velocity parameter (o"k), and the central phase-space density estimate (where a = \/3a K ). 
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Table 2. NCs parameters. 



Nuclear cluster 




M 


c 




r c 






Po 


tch 


O"0 


Po/^o 


NCi 


6.4 


x 10~ 2 


2.2 


6.3 


X 


io- 3 


0.12 


2.3 x 10 3 


0.16 


0.90 


3.1 x 10 3 


NC 2 


6.4 


x 10~ 2 


1.7 


5.0 


X 


io- 2 


0.21 


15 


0.38 


0.50 


120 



Note. — Parameters list for the NCs in the last configuration, expressed in galactic units. 
Reported are: the total mass (M) , the concentration coefficient (c) and the King radius (r c ) 
of the best-fit King model, the half- mass radius (rh), the central density (po), the half- mass 
crossing time (t c h), the central velocity dispersion (a ) and the central phase-space density 
estimate. 
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Table 3. Orbital initial conditions for the clusters. 



cluster model 


x 


yo 




(al) 


-0.4 


0.25 


-0.35 


(a2) 








(bl) 


0.3 


-0.35 


0.31 


(b2) 








(cl) 


0.375 


0.35 


-0.325 


(c2) 








(dl) 


-0.26 


-0.225 


0.425 


(d2) 









Note. — Initial conditions of the orbits of 
the 4 clusters. All started with zero initial 
velocity. 



